Multiple light scattering on the F = 1 — > F' = transition 
in a cold and high density 87 Rb vapor 
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We report an experimental study of near resonance light scattering on the F — 1 — > F' = com- 
ponent of the D2 line in atomic S7 Rb. Experiments are performed on spatially bi-Gaussian ultracold 
gas samples having peak densities ranging from about 5 • 10 12 — 5 • 10 13 atoms/cm 3 and for a range 
of resonance saturation parameters and detunings from atomic resonance. Time resolution of the 
scattered light intensity reveals dynamics of multiple light scattering, optical pumping, and satu- 
ration effects. The experimental results in steady-state are compared qualitatively with theoretical 
models of the light scattering process. The steady-state line shape of the excitation spectrum is in 
good qualitative agreement with these models. 

PACS numbers: 42.25.Dd, 42.50.Nn, 42.50.-P, 72.15.Rn, 37.10.Gh 



I. INTRODUCTION 

The interaction of light with dense atomic gases is a 
vigorous area of research in quantum optics [ll-[24| . Part 
of this interest stems from the interdisciplinary nature 
of the field, and with the large number of fundamentally 
important results and potential applications that have 
emerged. Recent efforts have ranged from basic studies of 
light localization in disordered systems [25:'- 33] , including 
atomic gases, to investigation of cooperative scattering 
[34l - l4lj ]. Areas of experimental and theoretical research 
with both fundamental motivations and possible applica- 
tions include searches for atomic physics based random 
lasing (42r!47! and quantum memories for quantum infor- 
mation and communications (48l - [50| . Single photon op- 
tical memories in dense atomic gases may potentially be 
derived from a number of approaches, including devel- 
opment of subradiant atomic-photonic modes [51] . and 
extensions of electromagnetically-induced-transparency 
(EIT) based approaches to novel two-photon optical 
schemes at higher densities 0, ©, [13, El IS SI SI • 

We have ongoing experimental and theoretical research 
efforts focused in part on developing quantum memories 
using either two-photon EIT based approaches on one 
hand [52T - [56l | . and on the possible formation of subradi- 
ant single photon modes on the other [27T - l29l | . To obtain 
formation of subradiant modes is a challenging experi- 
mental enterprise, and requires atomic densities ~ 10 14 
atoms/cm 3 , in order to achieve high orders of multiple 
light scattering in the sample. In addition, dynamical 
processes which either dephase the multiply scattered 
light or which lead to reduction in the scattering cross- 
section should be well understood. The current program 
is focused on nearly optically closed hyperfine transi- 
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tions associated with the D 2 transition in ultracold 87 Kb. 
These two transitions are the F = 2 — > F 1 = 3 main 
optical trapping component and the F = 1 — > F' = 
transition arising from the lower energy ground state hy- 
perfine component. We have reported elsewhere our ex- 
perimental research associated with light scattering on 
the F = 2 ->• F' = 3 transition [13 ■ 

In the present paper we report experimental and the- 
oretical investigation of the F = 1 — » F' = transi- 
tion in 87 Rb. Overall, our studies include examination 
of the roles of atomic density, optical saturation, and 
detuning of probe radiation from optical resonance on 
the F = 1 — » F' = transition. We will see that, for 
this transition, the debilitating effects of Zeeman optical 
pumping play an essential role in all aspects of the exper- 
iments. In the present paper, we concentrate primarily 
on steady state optical excitation and how the light scat- 
tering signals depend in that case on probe detuning and 
atomic density. In the following sections we first pro- 
vide some details of our experimental approach. This is 
followed by presentation of our experimental results and 
comparative discussion on the basis of theoretical models 
of the light scattering process. 



II. EXPERIMENTAL CONFIGURATION 

A schematic diagram of the experimental apparatus 
used in the measurements is shown in Fig. 1, while the 
particular optical transitions of interest associated with 
the atomic S7 Rb D2 line are shown in Fig. 2. As il- 
lustrated in Fig. 1, the central part of the experimen- 
tal apparatus is a magneto optical trap (MOT) which 
serves to form and confine cold s7 Rb atom samples. The 
MOT is a standard vapor-loaded trap formed in a vac- 
uum chamber with a base pressure ~ 10~ 9 Torr. The 
six MOT beams are derived from a single external cav- 
ity diode laser (ECDL) with the grating arranged in a 
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FIG. 1: Schematic diagram of the experimental apparatus. In 
the figure, PMT refers to an infrared sensitive photomultiplier 
tube and Amp refers to a fast preamplifier. MOT stands for 
magneto optical trap, while QUEST is an abbreviation for 
quasi electrostatic trap. 
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FIG. 2: Energy level diagrams illustrating the experimental 
scheme, (a) The MOT and repumping transitions used in 
the present experiment, (b) Probe excitation with linearly 
polarized light (z-direction) used in the present experiment. 
In part (b) of the figure, the Zeeman levels have been broken 
out to show the role of Zeeman optical pumping out of the F 
= 1, M = state to the F = 1, M = ± 1 states. 



Littrow configuration. The main master diode laser is 
frequency locked to a saturation absorption feature pro- 
duced in a room temperature Rb vapor cell. The master 
laser power is increased by injecting the master output 
into slave laser. The arrangement provides more than 
20 mW of trapping light in laser beams of cross sec- 
tional area ~ 2 cm 2 . The slave laser output is switched 
and spectrally shifted as required with an acousto op- 
tical modulator (AOM) to a frequency set at about 18 
MHz below the 87 Rb F = 2 -> F' = 3 trapping transi- 
tion. The repumper laser is also an ECDL of the same 
basic design as the main MOT laser, and is locked to the 
F = 1 — > F 1 = 2 hypcrfine transition. The repumper 
delivers a beam of maximum intensity ~ 0.6 mW/cm 2 
and is delivered along the same optical path as the main 
trapping laser beams. Switching of the repumper laser is 



also controlled with an AOM. 

In the experiments reported here, the cold atom sam- 
ple is initially produced in the higher energy F = 2 level. 
Direct absorption imaging measurements of the peak op- 
tical depth on the F = 2 — >• F' = 3 transition yielded, 
for this sample, b ~ 10 in a Gaussian radius of r Q ~ 
0.45 mm. However, the main sample production goal 
is to transfer a significant number of the trapped atoms 
to a carbon-dioxide laser (CO 2) based far off resonance 
optical dipole trap. The 100 W CO2 laser based trap 
operates at a wavelength of 10.6 fim and is deeply in the 
quasistatic trapping regime. This laser is focussed to a 
radial spot size of ~ 55 fim, and a corresponding Rayleigh 
range of zr ~ 750 \xm. The CO2 laser focal zone is over- 
lapped with the MOT trapping region, while application 
of the laser beam itself is controlled by a 40 MHz AOM. 
The atom sample formed in the MOT is compressed and 
loaded into the quasistatic dipole trap (QUEST) by de- 
tuning the MOT master laser 60 MHz to the low fre- 
quency side of the trapping transition, while simulta- 
neously lowering the repumper intensity over an order 
of magnitude. The resulting temporal dark spot MOT 
loads the atoms predominantly into the lower energy F 
= 1 hyperfine component. The result of this procedure is 
transfer of about 15 % of the MOT atoms to the QUEST. 
It is important to note that this transfer efficiency is de- 
termined after a QUEST holding period of about 1 sec- 
ond, during which the atomic sample naturally evolves 
towards thermal equilibrium (this happens through elas- 
tic collisions between the confined Rb atoms). Auxiliary 
measurements of the QUEST principal characteristics af- 
ter the 1 s hold period, by absorption imaging, paramet- 
ric resonance [581 ]. and the measured number of atoms 
transferred show a sample with peak density about 5 • 
10 13 atoms I cm 3 and a temperature of~ 65 fiK. The 1/e 
lifetime of the confined atoms is longer than 5 s, and is 
limited by background gas collisions. The residual mag- 
netic field in the sample area, when the MOT quadrupole 
field is switched off, is estimated to be less than a few mG. 

In the main experimental protocol, a probe beam 
tuned in the spectral vicinity of the F = 1 — > F' = 
nearly closed transition is directed towards the sample, 
and the resulting scattered light signals collected as il- 
lustrated schematically in Fig. 1. The probe laser is of 
the same design as the repumper laser, has a bandwidth 
~ 3 MHz, and is switched and directed by an acousto 
optical modulator towards the sample. Because of con- 
straints on the vacuum chamber geometry, the linearly 
polarized probe beam is directed (see Fig. 1) at an an- 
gle of approximately 30 degrees away from the fluores- 
cence collection direction. The probe beam is also di- 
rected downwards towards the sample at an angle of 22.5 
degrees (as shown in Fig. 1). Finally, the sample fluo- 
rescence is collected without regard to light polarization; 
as light scattering is dominated by the F = 1 — > F' = 
transition we expect the scattered light to be mainly un- 
polarized (very small contributions from the quite far off 
resonance F = 1 F' = 1,2 transitions are in general 
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A. On-resonance scattering 



TABLE I: QUEST parameters relating the peak transverse 
optical depth on the F = 1 — s> F' = transition to the 
maximum sample density and the Gaussian radii of the atomic 
cloud. 



polarized in the single scattering limit). 

To close this section, we point out that in some of the 
experiments reported here the atomic density was varied 
over a factor of about 10. This was accomplished by al- 
lowing for a period of ballistic expansion of the cloud after 
the QUEST was turned off. The atomic sample temper- 
ature is known (by ballistic expansion measurements), so 
this procedure allows the peak density or the peak optical 
depth to be determined. As the sample is well approxi- 
mated by a two-axis Gaussian atom distribution [59( , the 
two Gaussian radii and the peak atom density (or the 
total number of atoms in the sample), are sufficient to 
determine the two peak optical depths characterizing the 
sample. We summarize in Table 1 the peak transverse 
optical depth b t the peak atom density at the center of 
the sample n , the transverse Gaussian radius r a , and 
the longitudinal Gaussian radius z a . The optical depth 
refers here to that of the nearly closed F = 1 — > F' = 
hyperfine transition, which has a total resonance light 
scattering cross section of 3.23 x 1CU 10 cm 2 . 



III. RESULTS AND DISCUSSION 



In this section we present results associated with on- 
resonance light scattering, where on-resonance refers to 
the bare-atom F = 1 — » F' = hyperfine resonance fre- 
quency, and with variations around that frequency. We 
emphasize that the measurements are made at a fixed 
single angle with respect to the incident probe laser; the 
geometrical setup is described in the previous section. 
With this in mind, we first expect that the scattered 
light will be unpolarized, as the excited level has F' = 
0. However, because of the high optical depth of the 
sample, the scattered light intensity should show impor- 
tant angular dependence [54| [60] ■ For instance, light 
scattered in the near forward direction shows a minimum 
in the spectral variations near the resonance line center. 
This can effect the time evolution of the scattered light 
signals as well 54j . On the other hand, light scattered 
in the backwards direction shows an enhancement due 
to the coherent backscattering effect. This spatially cone 
shaped feature depends critically on probe laser inten- 
sities near saturation and also on the spectral detuning 
from atomic resonance. 
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FIG. 3: Time evolution of the probe scattered light intensity 
on the F — 1 — > F' = hyperfine transition. The atomic 
density is at its peak level for these measurements. Signals 
are shown for several different probe laser intensities. The role 
of optical pumping is apparent at the larger probe intensities. 
Probe detuning from resonance A = 0. 
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FIG. 4: Variation of the peak and steady state signal levels 
as a function of the probe laser intensity. Atomic density is 
at is maximum value n for these data. Probe detuning from 
resonance A = 0. 

We start by presenting in Fig. 3 the results of typical 
measurements of the time dependence of light scattered 
from the atomic sample. These measurements are made 
for a maximum density sample, which in this case has a 
transverse optical depth of bt ~ 40, and a peak atom den- 
sity at the center of the sample of 5.0 • 10 13 atoms/cm 3 . 
These measurements are made for a range of probe laser 
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intensities, which for the most part are well below the 
saturation intensity I sat ~ 14 mW/cm 2 for linearly po- 
larized excitation of the F = 1 — >• F' = transition. 
We see for higher probe laser intensities a quite rapid 
buildup to a peak intensity followed by a decay of about 
1 /Lts to a steady state. The build up time is less than 
about 100 ns, and is limited by the turn-on time of the 
switching AOM used to control application of the probe 
laser to the sample. We also note that as the probe laser 
intensity is reduced, the transient peak is reduced, and 
for the lowest probe laser intensities the scattering signal 
smoothly rises to a steady value. The overall behavior is 
summarized in Fig. 4, where the nearly linear growth of 
the peak in the scattering signal, and the leveling off of 
the steady state signal is shown. This general behavior of 
the transient response and the variations of the response 
with probe laser intensity was found to be qualitatively 
the same for a range of densities from 0.1 n D up to the 
peak density n Q (see Table 1). 
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FIG. 5: Variation of the total scattered signal intensity due 
to systematic changes in the atomic density. Probe detuning 
from resonance A = 0. 

As the maximum laser intensity used in these experi- 
ments is well below the optical saturation intensity, we 
can understand the transient behavior to result from Zee- 
man optical pumping in the ground level due to elas- 
tic Raman transitions to the F — 1,M — ±1 states. 
This optical pumping does not completely deplete the 
F = 1, M = state because multiple elastic light scatter- 
ing allows for a build up of radiation within the atomic 
sample. This gives rise to optical pumping by the dif- 
fuse light in the sample, which competes with direct Zee- 
man optical pumping by the probe beam and results in 
a nonzero steady state level for the scattered light signal 
as seen in Fig. 3 and Fig. 4. This process as consid- 
ered for a wide range of densities and laser intensities 
will be treated in a later report. In the present paper we 
are concerned with the steady state signals for the case 
where the probe intensity is so weak, and the number 
of scattered photons so few, that there is no mesoscopic 



rearrangement of populations in the F = 1 level during 
realization of a single sample. The main results of this 
paper, as presented in this and the following section, are 
recorded under those conditions. 

We present in Fig. 5 measurements of the variation 
of the total scattered light intensity from the F = 1 — > 
F' = hyperfme transition as a function of systematic 
changes in the atomic density. These measurements were 
made at very low probe laser intensity, under conditions 
selected so that there was no perceptible optical pump- 
ing, as in the lowest intensity results shown in Fig. 3 and 
Fig. 4. In Fig. 5 we see that as the density decreases, 
the overall intensity of the scattered light increases. A 
similar behavior [53, [6(| has been observed for measure- 
ments on the F = 2 — > F' = 3 hypcrfinc transition of 
87 Rb. The effect is due to the collective nature of near 
resonance light scattering from a high density and cold 
atomic gas. Because the optical depth is so large for the 
highest density, light scattering occurs primarily from the 
outer surface of the sample. Relatively few atoms then 
contribute to the sample. For lower density the light can 
penetrate more deeply and a relatively larger number of 
atoms participate, leading to a larger signal. 



B. Variations with spectral detuning from 
resonance 

We now consider the variations of the scattered light 
signals as a function of detuning of the probe laser fre- 
quency from bare atomic resonance, where A = 0. These 
data are recorded under conditions of a very weak probe 
laser, so that Zeeman optical pumping is negligible. For 
the steady state regime (as in the lowest probe intensity 
results of Fig. 3), this dependence is shown in Fig. 6. 
There we see that the spectral variations show a clear 
resonance behavior in a range of 3 to 4 7 about A = 0, 
where 7 is the ~ 6 MHz natural width of the atomic 
resonance. The solid line is a Lorentzian line profile, 
which well describes the resonance line shape, and yields 
a full width at half maximum of slightly larger than 10 
MHz. A full-width greater than the natural width is 
expected, due to absorption broadening and due to the 
dipole dipole interaction of the atoms under high density 
conditions. Selection of the Lorentzian line shape is arbi- 
trary, but gives a decent fit to the data and a consistent 
way of estimating the full width under different experi- 
mental conditions. Within the spread of the data, a fit 
using a Gaussian line profile yields essentially the same 
full width at half maximum. 

To further analyze the experimental data, we have ex- 
tracted the full width at half maximum of the spectral 
profile as a function of time [57]. With reference to the 
lowest intensity results of Fig. 3, this is equivalent to 
making a slice on the time axis of 100 ns and recording 
the spectral profile as a function of detuning in this time 
window. Repeating this procedure for the full sequence 
of data over the 2.5 /zs range of the probe excitation and 
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FIG. 6: Near-resonance spectral variation of the total scat- 
tered light intensity in the spectral vicinity of the F = 1 — > 
F' — hyperfine transition. 
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FIG. 7: Time evolution of the probe scattered light intensity 
on the F — 1 — > F = hyperfine transition. The atomic 
density is at its peak level for these measurements. 



decay signals yields the excitation spectrum shown in Fig. 
7. There we see that the line width is very large for short 
times, rapidly decays to its steady state level (as in Fig. 
6), and then sharply decreases again after the probe pulse 
is extinguished. For the short-time turn on of the excita- 
tion pulse single scattering dominates, and the line width 
can be estimated to be A Q = j/2y/b a — 1, where b is the 
on resonance transverse optical depth through the center 
of the cloud. For the experiments reported here b D = 40, 
giving a full width of A Q = — 1 ~ 37 MHz, in very 

good agreement with the short time value of Fig. 6. 

Upon turn off of the probe pulse the spectral width of 
the excitation spectrum decreases from its steady state 
value of about 10 MHz to a value on the order of the 
natural width. As in the spectral response upon turn-on 
of the prober pulse, the line shape here is also well fit by 
a Lorentzian form. However, we point out that the spec- 
tral width of the probe laser itself is about 1 MHz, and 
has a measured Gaussian power spectrum. This means 
that the widths determined by these measurements are 
slightly larger than that determined by the physical pro- 
cesses involved. In fact, the longest lived mode [13] for 
these samples corresponds to the so called Holstein mode, 
or the longest lived diffusive mode for the sample un- 
der study; the lifetime of this mode can be significantly 
smaller than the natural width of the transition. 



IV. THEORY 

In this section we present theoretical calculations for 
comparison with the steady state line shape of the exci- 
tation spectrum. In general, microscopic calculations for 
the cooperative scattering process for a macroscopic col- 
lection of atoms with a degenerate ground state is rather 



difficult because of the rapidly rising number of equa- 
tions to be solved deNd^ -1 , where d e is the degeneracy 
of the atomic excited state and d g is the degeneracy of the 
ground state. Here N is the number of atoms considered 
in the calculations. This makes practically impossible 
exact microscopic calculations for the scattering process 
and for its time-dependent fluorescence dynamics. In this 
section we present the results of our calculations of the 
scattering spectrum in the steady state regime, which are 
based on two complementary models of the self-consistent 
and approximate microscopic approaches, see Ref. [27j . 
and discuss the results of our numerical simulations in 
the context of the data presented in the experimental 
part of the paper. 

The self-consistent approach accepts the description 
of the scattering problem on the level of the macroscopic 
(i.e. mesoscopically averaged) Maxwell theory, where the 
dielectric susceptibility obeys an equation resulting from 
the self-consistent dynamics of the driving (probe) field 
and of the atomic dipoles. The crucial assumption is 
that groups of closely located atomic dipoles respond to 
the field cooperatively such that the longitudinal dipole- 
dipole interaction can be incorporated into the Lorentz- 
Lorenz local-field correction. For a low intensity of the 
driving field (where optical pumping can be neglected) 
the atoms equally populate the Zeeman sublevels and the 
medium is isotropic and can be parameterized by a sin- 
gle dielectric constant. If the atoms are homogeneously 
distributed in a spherical volume, then it is possible to 
express the scattering cross section by the standard solu- 
tion of the Debye-Mie problem. In Fig. [3] we reproduce 
the spectral dependencies of the total cross section calcu- 
lated for such a spherically symmetrical atomic system. 
The chosen radius of the sphere 30 fim is roughly scaled 
by the spot size where the probe beam crosses the atomic 
sample. We use the densities essentially less than in the 
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FIG. 8: The scattering cross section calculated in the Debye- 
Mie model for an atomic sample with radius 30 /im. Atoms fill 
a sphere homogenously with the density varying from 0.005 
to 0.025. 



peak point of the cloud to emphasize the importance of 
the contribution from the tail area of the Gaussian dis- 
tribution in the experimental observation of the scatter- 
ing process. The plotted dependencies demonstrate the 
smoothed spectral profile with a bandwidth qualitatively 
in agreement with the experimental data of Fig. [Sj Note 
that there is a slight asymmetry in the spectral behavior 
related with the density effects which is a consequence 
of the spectral asymmetry of the dielectric permittiv- 
ity. Because of the essential differences in geometries of 
this model and the experiment, which is a prolate dou- 
ble Gaussian distribution in experiment and spherically 
homogeneous in theory, we cannot make direct quantita- 
tive comparison of experimental and theoretical results 
but we can point out at least the qualitative agreement 
between both the data sets. We finally point out that 
the near resonance theoretical behavior is structured as 
a competition of longitudinal red shift and cooperative 
radiation blue shift but this effect is very likely masked 
by the uncertainty in experimental data. 

We note that the time dynamics associated with the 
fluorescence decay can be also modeled in the framework 
of the self-consistent calculation scheme. However, this 
would require the complete Monte- Carlo simulations of 
the process with the Green's function formalism, see Ref. 
[2^ | . The scattered energy of the light scattered from a 
dense atomic cloud in the steady state regime emerges 
mainly from its surface area and the transient deviations 
for the widths of the fluorescence spectra reproduced in 
Fig. [7] could be verified by such a Monte- Carlo simulated 
and spectrally sensitive self-consistent calculations. 

As we commented above, the microscopic calculations 
can be only approximately done for the considered tran- 
sition. Here we reproduce the results of our calculations 
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FIG. 9: The scattering cross section calculated microscopi- 
cally for an ensemble consisting of one hundred atoms. Atoms 
have a spherically symmetrical Gaussian distribution with the 
peak density varying from noA 3 = 0.001 to 0.2 and populating 
only one Zeeman sublevel. 



based on the general formalism of the quantum scattering 
theory previously developed in Ref. [27| and present the 
data for the case of a sample consisting of one hundred 
atoms. For the typical experimentally attained atomic 
densities, which are less than rtoA 3 < 0.1 in its peak value, 
in order to resolve the quasi-energy structure, associated 
with the longitudinal dipole-dipole interaction, it would 
be enough to take into consideration only a few nearest 
neighbors in the vicinity of each selected scatterer. This 
scheme were tested by us on a small number of atoms 
and it demonstrated reliable convergence, which reduces 
the number of equations to be solved to deNd™ -1 , where 
n is the number of the neighboring atoms retained in the 
calculation. 

In Fig. [5]we show the results of our calculations, which 
were done for n = 4 and for all the atoms populating 
one Zeeman sublevel. Let us point out here again that 
because of high degeneracy of the ground state our results 
cannot be statistically averaged over all the populated 
initial states. That means we cannot directly compare 
them on a one-to-one basis either with experiment or 
with the self-consistently calculated data. Instead, we 
rely on qualitative comparison among the results. 

First, as was explained in Ref. [27j, because of the 
presence of elastic Raman scattering channels there is no 
visible signature of either super- or sub-radiant Dicke- 
type exciton modes for the F = 1 — > F' = transition. 
In contrast with F = — > F' = 1 transition (where 
such states are normally predicted and discussed) in the 
F = 1 — > F' = case all the resolvent poles are described 
as the resonances with a typical line width of about 7. 
Thus the presence of the cooperative longitudinal and ra- 
diative interactions manifest themselves in the smoothed 
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variation of the scattering spectrum associated with dis- 
order and with the configuration dependence. The slight 
modulation of the experimental spectrum, which is sug- 
gested in the experimental data shown in Fig. [5] can 
be not only experimental uncertainty but can be also 
manifestation of the cooperativity and disorder effects. 
Second, we can point out that because of competition 
between near and far field interactions the microscopi- 
cally calculated spectra have an asymmetry in respect 
to the atomic resonance line. The latter observation is 
in qualitative agreement with the predictions of the self- 
consistent model, which were pointed out above. 

V. CONCLUSIONS 

We have presented detailed experimental and theoret- 
ical results associated with scattering of light from an 
ultracold and highly dense gas of 87 Rb atoms. For radia- 
tion tuned in the spectral vicinity of the F = 1 — > F' = 
hyperfine component of the Di line, we have studied the 
atomic density, probe laser detuning and probe laser in- 
tensity dependence of the scattered light intensity. The 
measured time dependence and steady state responses in- 
dicate that the light dynamics is strongly effected by Zee- 
man optical pumping in the lower energy F = 1 hyperfine 
component. However, for very low probe laser intensities, 



Zeeman optical pumping can be made negligible; this has 
allowed for study of the steady state regime in the ab- 
sence of optical pumping. Good results are obtained in 
comparison between the experimental results and theo- 
retical calculations. The strong role of elastic Raman 
transitions also strongly suggests that this transition is 
not particularly suitable for searches for Anderson light 
localization in three dimensions. 
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